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1.0 INTRODUCTION 


Under NASA Contract NAS8-33541 , a prototype of a fog dispersal unit 
was built and preliminary performance tests carried out. A general 
description of the system and discussion of the results are presented in 
the final report from that contract (Frost 1982). A major conclusion 
by Frost (1982) was that a number of significant design modifications 
need to be investigated. These modifications have been made and follow- 
on testing has been carried out. Beneficial experience toward the 
design of both a second-generation operational -type unit and of an 
advanced research-type unit has been obtained from these tests. The 
major modifications were made to the nozzle and to the flow control 
system. Moreover, a new system of instrumentation to record the current 
output and to measure the voltage drop at the corona was developed. 

This report describes the modifications to and subsequent testing 
of the charged particle generator. Section 2.0 describes the hardware 
modifications to the experimental apparatus. Section 3.0 presents the 
results and compares these with previous results reported by Frost 
(1982). The results are interpreted in terms of improvements to the 
performance of the system. Section 4.0 presents the conclusions rela- 
tive to the modified prototype unit. 
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2.0 EXPERIMENTAL APPARATUS 


Modifications to the previous fog dispersal unit under this contract 
effort are reported in this section. Each modification is described and 
compared with the previous system arrangement. The description of the 
experimental apparatus was divided into five sections in the last 
report. These sections were: (1) air supply circuit, (2) nozzle, 

(3) water flow circuit, (4) electrical circuit, and (5) current output 
measurement circuit. The modification to each of these categories is 
described in the following subsections. 

2.1 Air Supply Circuit 

The general arrangement of the experimental apparatus is illus- 
trated in Figure 2.1. The previous arrangement used the compressor to 
supply the airflow through the nozzle and to act as a high-pressure 
source to drive the water flow. However, this arrangement resulted in 
poor control of the water flow rate. Consequently, a separate positive 
displacement pump was added to the system to pump the water, and the 
compressor was used only to supply airflow through the nozzle. Other 
than this modification, the air supply circuit remains the same. 

2.2 Nozzle 

A new converging/diverging nozzle was cast out of Plexiglas. This 
nozzle was cast with an attractor at the throat and was designed for an 
exit Mach number of 1.35. The nozzle used on the original prototype 
system was also designed for a Mach number of 1.35; however, due to the 
necessity of using a thicker corona needle than originally planned, the 
nozzle did not perform as designed. The geometry of the new nozzle was 
designed to account for this blockage effect. In actual fact, the 
geometry of the new nozzle is not different in dimensions from that of 
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Figure 2.1 Schematic of charged particle generator 




the previous nozzle because the attractor is now at the throat and the 
needle does not protrude into the throat region, and hence does not 
reduce the actual throat area. 

The cross sectional area of the throat controls the Mach number at 
the exit for a given pressure difference across the nozzle. Again, 
however, difficulties were encountered with this new nozzle. During 
casting of the nozzle, a good seal between the wire running to the 
attractor and the Plexiglas was not achieved. Consequently, water and 
air leaked through this region. Arcing thus occurred at high-voltage 
settings and eventually burned out the wire. Consequently, the nozzle 
was not useful for further testing, and all results reported herein are 
for tests utilizing the older nozzle. 

2.3 Water Flow Circuit 

Liquid water is supplied from a 12-gallon tank. The tank was 
originally pressurized with bleed airflow from the primary airflow used 
to supply the nozzle. With this arrangement, extreme difficulty in 
controlling the liquid water flow rate was encountered. Therefore, a 
FMI* laboratory pump RP-SY-1 was installed in the water supply circuit 
and the pressure bleed airflow arrangement disconnected. The water is 
now pumped from the water tank through a rotameter to the atomizer valve. 
Figure 2.2 is a photograph of the water pump. The RP-SY-1 laboratory 
pump can control the flow rate from 1.9 cc/min to 19.4 cc/min and can 
attain head pressure up to 60 psi. 

2.4 Electrical Circuit 

A 20,000 volt power supply with a variable rheostat control circuit, 
as described in Frost (1982), was used for all tests carried out in this 
study. The voltage drop across the attractor and needle was measured 


*Trade name of FMI Fluid Metering, Inc., 29 Orchard Street, 
Oyster Bay, NY 11771. 
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Figure 2.2 FMI RP-SY-1 water pump. 

with a B&K** 3-1/2-digit digital multimeter connected with a B&K PR-28 
high-voltage probe (see Appendix A). The probe was placed in the 
circuit, as shown in Figure 2.3. Using the PR-28 high-voltage probe 
coupled with the digital multi-range voltmeter, voltages up to 40 kV DC 
can be measured to ±3 percent accuracy. The high-voltage probe proved 
very successful in measuring both the voltage drop across the corona and 
the voltage along the centerline of the jet. 

2.5 Current Output Measurement Circuit 

Two output current measuring systems were used for the present 
tests. One was the previously used grounded needle arrangement 

**Trade name of Dynascan Corporation, 6460 West Cortland Street, 
Chicago, IL 60635. 
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Figure 2.3 Circuit for measuring the voltage drop across the gap 
between the corona attractor and needle. 


(Frost 1982). The other was a current collection system, based on 
discussions with Lothar Ruhnke of the Naval Research Laboratory. This 
system consisted of a conductive cylinder placed over the jet to collect 
all charges exiting the nozzle. The cylinder was grounded through an 
ammeter, which measured the current in the jet. This cylinder was 6 
inches in diameter, 2 feet long, and covered on the top with an aluminum 
mesh. The configuration of this collection system is shown in Figure 
2.4. The following reported measurements of output current have been 
made with both the new collection system (labeled pipe measurement) and 
the previously used needle probe system (labeled needle measurement). 
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Figure 2.4 Charged droplet collection pipe. 


7 




3.0 RESULTS AND DISCUSSION 


Results of the tests carried out are discussed in this section. 

The output current with variation of the controllable parameters is 
the primary variable measured. Comparison with previous results is 
also given and the effects of modification to the previous system 
described. In addition, measured velocity and voltage distributions 
along the jet not made during the previous study are discussed. 

3.1 Current Output 

As in the previous report (Frost 1982), the controllable param- 
eters include: (1) the liquid water flow rate, (2) the corona voltage, 

and (3) the nozzle stagnation pressure. Figure 3.1 schematically illus- 
trates the measurement arrangement. The current output measurements 
include both those made with the pipe system and those made with the 
needle probe. Comparison is made with previous results. 

3.1.1 Liquid Water Flow Rate 

The variation of the electrical current of the jet with liquid 
water flow rate is shown in Figure 3.2. In this figure, the current 
output is scaled by the stagnation pressure and the relative humidity, 
RH, of the atmosphere. The output current was found to scale with RH 
empirically. This effect is believed to be associated with the depen- 
dence of the electrical conductivity of the atmospheric air on RH. 
Conductivity will increase with increasing RH allowing more current to 
leak to ground in preference to flowing through the measuring probe. 
Thus, under very dry conditions, a higher jet current measurement is 
anticipated. 

Results for both the needle measurement and the pipe measurement 
are shown on the figure for the same test. These measurements were made 
by changing probes without shutting down the system or changing the 
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Figure 3.1 Schematic of output current measurement with pipe and needle probe. 


Liquid Jet Current (I»RH/P-j) 



Figure 3.2 Comparison of the variation of scaled current output versus liquid water 
flow rate with two different measurement techniques. 
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control settings. The bottom of the pipe and the needle were placed 1/2 
inch above the plane of the nozzle exit. Both systems measured essen- 
tially the same values. Any difference is attributed to experiment 
variability. 

The optimum current output occurs at a liquid water flow rate in 
the range of 4 to 7 cc/min. Figure 3.3 compares the recent results with 
previous measurements. Both recent and previous results agree reasonably 
well. The scatter in the data is again attributed to experiment vari- 
ability; particularly in view of the difficulty associated with main- 
taining a stable water flow rate in previous tests. 

3.1.2 Corona Voltage 

The circuit for measuring the corona voltage is shown in Figure 
2.3. Figure 3.4 shows the measured variation of jet current output with 
the corona voltage. The output current is roughly proportional to the 
corona voltage over the range 4 to 12 kV. Generally, the voltage drop 
across the gap between attractor and the needle was not stable. This 
was due to the natural unsteady nature of a corona discharge (Loeb 1965). 
Also, the turbulent flow in the nozzle vibrates the needle and changes 
the electrical field of the corona, causing a corresponding change in 
the output current. 

3.1.3 Vertical Variation of Current in Air Je t 

The vertical variation of electrical current in the air jet was 
measured by moving the needle probe and pipe vertically. Figure 3.5 
shows the results of both the needle probe and pipe measurements. 

Again, there is no apparent difference. Figure 3.6 shows the comparison 
of the recent results and previous results. The recent results appear 
to be roughly 20 percent higher than the previous results measured with 
the needle. No effort to scale these results has been made. Note, 
however, the scaling in Figure 3.3 tends to collapse the data onto one 
curve. Other differences may be due to experimental variability for 
each run caused by small changes in the position of the needle relative 
to the attractor, for example, and by variation in cleanness of the air 
filter, etc. Variation of current with height results from charges 
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Figure 3.3 Comparison of the scaled current output from recent results with previously 
measured values. (H = height above the exit plane of the nozzle.) 



Figure 3.4 Variation in jet current output with corona voltage. 


migrating out of the jet due to turbulence and their own mobility in an 
electric field. 


3.2 Voltage Potential Distribution 

The variation of voltage potential along the air jet was measured 
by adjusting the needle probe vertically and connecting the needle probe 
in series with the high-voltage probe and digital voltmeter, as shown in 
Figure 3.7. Figure 3.8 shows the measured vertical voltage distribution 
in the air jet for two different values of the corona supply voltage. 

The trend of both sets of data are the same. The voltage potential 
increased from the exit of the nozzle and reached a maximum at 2 to 3 
inches above the exit and then decreased gradually. 

The peak in the voltage distribution in the vertical direction is 
associated with the point where the kinetic energy of the droplets 
issuing from the nozzle is converted to electrical potential energy. 
Above the height at which the peak voltage occurs, the charged droplets 
are displaced upwards due to their mutual repulsion. 
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Figure 3.5 Comparison of the jet current versus height above exit plane of the nozzle 
for the two different measuring techniques. 
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Figure 3.6 Comparison of the measured current output versus height above the exit 
plane of the nozzle with previous measurements (Frost 1982). 
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Figure 3.7 Illustration of jet voltage measurement. 
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Figure 3.8 Vertical variation of voltage along the charged air jet (water flow 
rate m = 6 cc/min) . 



Ruhnke (1970) carried out a computer analysis where he solved both 
the Poisson electric field equation and the continuity equation for 
current. Figure 3.9 is a plot of his numerical results. The data mea- 
sured in these tests are compared with the theoretical model by scaling 
the voltage with the maximum voltage and the height with the height at 
which the maximum voltage occurs, z . Figure 3.10 compares the measured 
data with the numerical results. The experimental data are found to 
bracket the analytical results very well. 

3.3 Velocity Distribution in Air Jet 

A pitot tube was used to measure the jet velocity in both the 
vertical and radial directions. The ALNOR^ 6000-P velocimeter, which is 
generally used for measuring low-speed channel flow velocities, was used 
for this purpose. The range of this pitot tube is only 0 to 10,000 fpm 
(0 to 50.8 m/s). Thus, the meter read off-scale when measuring the jet 
flow velocity from the exit of the nozzle to approximately 10 inches above 
the plane of the nozzle. Above this height the meter provided meaning- 
ful results. Figure 3.11 shows the measured vertical velocity profiles 
at different heights above the nozzle exit plane. While the measured 
data is rather rough because of the crosswind effect, the trends of the 
distributions are clearly indicated. As expected, the velocity profile 
is sharper and narrower at lower levels and flatter and wider at higher 
levels. Figure 3.12 shows the variation of the centerline velocity of 
the jet with height. Again, the data is measured only at heights above 
10 inches from the jet exit plane. 

Measurements of the jet velocity were compared with theoretical 
models and other empirical models reported in Rajaratnam (1976). 

Rajaratnam predicts that the velocity along the centerline of the jet 
scaled with the velocity at the jet exit varies inversely with the 
distance along the jet scaled with the jet diameter. However, the 
axial decay of the velocity for circular jets in cross flow, as predicted 


tTrade name of Alnor Instrument Company, 7301 North Caldwell 
Avenue, Niles, IL 60648. 
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Figure 3.9 Computed voltage potential distribution in a space charae 
jet (Ruhnke 1979). y 



Figure 3.10 Comparison of scaled measured data and numerical predictions 
(Ruhnke 1979). (z m = height above nozzle exit plane where 
voltage is maximum.) 
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by Patrick (1967), is illustrated in the insert of Figure 3.13. The 
measured data shown in Figure 3.13, similarly scaled, clearly illus- 
trates the same behavior as that of the circular jet in cross flow. 

For working purposes, an empirical curve fit of the data was com- 
puted as shown in Figure 3.14. The scaled velocity is found to corre- 

-1 5 

late roughly with z/d Q ' . The results of comparing the measured jet 
velocity with previous experimental correlations is believed to indicate 
that the velocity measurements, although relatively crude, provide mean- 
ingful results. 


3.4 New Nozzle Evaluation 

The new nozzle was cast with the attractor at the throat because it 
was anticipated that condensation occurred at that location. Unfortu- 
nately, there were some defects in the new nozzle. Without the needle 
blockage at the throat, the throat area was now too large for the 15 
SCFM compressor to provide sufficient flow to achieve choked conditions. 
Thus, Mach 1.00 at the throat did not occur, greatly reducing any drop- 
let condensation which might occur there (Wegener and Mack 1958). The 
attractor and electrical lead were not properly sealed to the casting 
with the result that there was flow leakage through the gap between the 
electrical lead and the Plexiglas. At high corona voltages, severe 
arcing occurred at this spot and resulted in burning of the wire insula- 
tion, compounding the problem. Eventually the nozzle could no longer be 
used. 

The new nozzle was tested but did not perform as expected because 
of the above-mentioned difficulties. The maximum stagnation pressure 
was 24 psig with the pressure regulator completely opened. Higher pres- 
sure could not be reached because the compressor could not supply a higher 
flow rate. Under these airflow conditions, the output current was only 
8 ya when the rheostat was set to approximately 8,000 volts. 

The new nozzle thus provided little (some quantitative effects were 
observed but the evidence is weak and they will not be discussed at this 
time) if any useful data. Consequently, all reported results are for 
tests carried out with the old nozzle. 
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Figure 3.13 Measured axial decay of scaled velocity (insert from 
Patrick 1967). (e; = curved distance along the jet 

centerline. ) 
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Figure 3.14 Curve fit on measured centerline jet velocity versus height 
above nozzle exit plane. 
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To correct the above-noted deficiencies in the nozzle, it is 
necessary to machine a new centerpiece for the mold used to cast the 
nozzle. The machining of the centerpiece is quite intricate and will be 
relatively expensive. To improve the seal between the attractor wire 
and the casting, it will be necessary to pour a new casting. This is 
not expensive, but time-consuming. Casting of the nozzle is not routine 
but requires some artistry. It is believed, however, that this will not 
create a problem in mass production of the electric charged particle 
generators since, for this operation, a production-type mold and casting 
procedure will be developed. 
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4.0 CONCLUSIONS 


The results of the experiment carried out in this study indicate 
that a charged particle generator, which will routinely generate an 
output current of 20 ya or greater, can be constructed and will operate 
consistently. Conclusions from the experiment relative to the proposed 
nozzle modifications were not reached because of defects in the nozzle 
which could not be overcome within the time frame of the study. The 
addition of a positive displacement pump to supply the liquid water to 
the atomizer valve proved highly successful. With this new arrangement, 
the charged particle generator, using the old nozzle, ran continuously 
without adjustments to the controls. Other conclusions reached during 
the study are as follows. 

1. Measurements of the jet output current using both the 
cylindrical pipe collector and the needle probe arrange- 
ment showed identical results within the variability of 
the experiment. It is concluded, therefore, that all 
previous measurements using a needle probe are meaningful 
and that the measured variation of electric current output 
with height along the centerline of the jet and as a func- 
tion of the liquid water supply are true measures of the 
charged particle generator output capabilities. 

Additionally, measurements of the voltage using a special 
high-voltage probe show the experimental data to behave 
in accordance with theoretical analyses based on numerical 
solutions of a simple turbulent jet flow model with the 
Poisson field equation and continuity of current. 

2. Measured velocities along the centerline of the jet com- 
pare well with other reported experimental results 
obtained in the laboratory and from theoretical analyses. 

The empirical curve fit of the centerline velocity provides 
a useful working equation for analyzing the ability of the 
jet to disburse charged particles to given height as 
described in Appendix B. 

3. Measurements of the corona voltage drop versus the output 
current from the jet show a roughly linear relationship 
from 4 to 12 kV. 
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Based on the measured results, an estimate of the size of the 
charged water droplets can be made. This analysis is given in Appendix 
C. The results indicate that the droplet size lies in the range of 0.8 
to 0.4 y. 

The charge to mass ratio computed for the prototype generator, 
i.e., output current divided by liquid water flow rate, gives an esti- 
mated result of 0*2 to 0.3 couls/kg. Lapple (1972) indicates that the 
range of charge to mass ratio achieved with a needle/attractor system as 
employed in this study is from 0.01 to 5 couls/kg. This observation is 
consistent with the measured results. 
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APPENDIX A 


SPECIFICATIONS OF HIGH-VOLTAGE PROBE 


The B & K-Precision high-voltage probe is a self-contained range multiplier 
designed to extend the voltage range of any constant input impedance type 
voltmeter. It is particularly suited for use with the following B & K-Precision 
instruments: 


MODEL NO. 


277 

282 

2800 

2830 

280 

2S3 

2810 


281 

290 

2815 



The outside case is made of high-impact ABS type plastic, which insures 
excellent insulation, impact resistance, and light weight. 


SPECIFICATIONS 

Voltage Range: 


Impedance: 
Accuracy: 
Calibration: 
Frequency Range: 


Usable to 40 kV DC. 

Multiplies instrument range by 1000. 

(E.g., 50V becomes 50,000V). 

May be used with any instrument 
having a constant input impedance 
of 10 or 15 megohms. 

600 megohms nominal. 

±3%, plus instrument accuracy. 
Internal; factory calibrated at 25 kV. 
DC to 80 Hz. 


IMPORTANCE OF MEASURING THE HIGH VOLTAGE 


All color TV sets require a specified value of high voltage If the high voltage is 
incorrect, the set may not function properly. Before any color set can be accurately 
converged, the high voltage must be checked and/or adjusted to the manufacturer’s 
specifications. 

A second reason for checking the high voltage arises from the possibility of 
X-ray radiation from certain sets, particularly older models. The possibility of this 
hazard is increased if the high voltage of the receiver is excessive. Always check the 
manufacturer’s instructions and specifications prior to performing high-voltage 
measurements. 
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APPENDIX B 


PREDICTION OF CHARGED PARTICLE TRAJECTORY 


The motion of a charged particle in a jet is governed by the drag 
due to the jet flow, the electrical field, and the gravitational field. 
Based on the assumption that the temperature, density, gravity, and 
number of charges on the particle are constant, the equations of motion 
in Cartesian coordinates gives: 


mx = 


my = 



a A c D (x -W x )|V a 
a AC D<* - V! V a 


+ F 


ex 


+ F 


ey 


(B.l) 


m'£= - ^ P a AC D (z - W z )iV a [ ♦ F ez 


- mg 


where m is the mass of the particle with cross sectional area A; x, y, 
z are the particle accelerations in the x, y, and z directions, respec- 
tively; p a is the air density; C^ is the drag coefficient defined as a 

function of Reynolds number; W , W , and W are the three instantaneous 

x y z 

wind velocity components; |V | is the magnitude of the particle velocity 
relative to the air; g is the acceleration of gravity; and F , F , and 

GX Gy 

F ez are the components of the electrical force. 

The three instantaneous wind velocity components W , W , and W , 

x y z 

which appear in Equation B.l, result from the air jet; this assumes 
no meteorological wind field. A free jet model from Rajaratnam (1976) 
is used in the model. Figure B.l is a sketch of a circular turbulent 
jet. The axial velocity of the jet in the fully developed region is 
given by 


u(r,z) 


m 


1 + 0.125 ? 2 


(B.2) 
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where % = ar/z and a - 18.5 and the radial velocity is given by 


v(r,z) 


U m^ z ^ x - 0.125 g 3 
0 2(1 + 0.125 C 2 ) 2 


(B.3) 


The variation of the centerline velocity, U m (z), is shown in Figure 3.14, 
page 22. The curve fit is given by 


U (z) = 
rrr ' 


1 + 0.014 


' z .1.5 


(B .4) 


where U is the nozzle exit velocity of the air jet. d^ is the diameter 
o o 

of the nozzle exit. Figure B.2 illustrates the axial velocity given by 
Equation B.4 at different heights. 

In addition to the jet velocity field, the electrical field must be 

specified. The three components of the electrical field, F ov , F , and 

ex ey 

^ez’ in Equation B.l can be expressed as 


F_„ = 1.6 x 10' 9 nE 
ex x 

F ey = K6 x 10' 9 nE y (B.5) 

F_ = 1.6 x 10' 9 nE 
ez z 

where n is the number of charges on each particle and E x , E y , and E z are 
the three components of the electrical field vector. From Ruhnke's (1970) 
analysis, the voltage potential field is as shown in Figure 3.9, page 19. 
The comparison between measured data and Ruhnke's (1970) calculation is 
discussed on page 18 and is shown in Figure 3.10, page 19. 

The electrical field strength along the center of the jet is plotted 
in Figure B.3. The data in Figure B.3 are interpolated from the results 
of Ruhnke's (1979) numerical analysis (see Figure 3.9). In Figure B.3 
the electrical field strength is normalized by E(0), which is the strength 

32 






i.Or 


See Eq. B.6b 



Figure B.3 Electrical field strength along centerline of vertical jet 
(z m is the vertical distance from nozzle exit to where 
vol tage is maximum) . 


at the nozzle exit. Also, the height is normalized by z m , which is the 
vertical distance from the nozzle exit to the point where the magnitude 
of voltage potential is a maximum. 

As Figure B.3 illustrates, the electrical field along the center- 
line of the jet can be divided into three regions. Curve fits of the 
electrical field for each region are 


E(z) 


when z < 


10 + (z/d Q ) 2 


- z m 


z - z_ 


O- 1 E(0 > 0JT" he " < 

m 


9 - 4. E(0) . when 1 .6 


1 + z/z m - 0.6 
m 


'm 


z < 1.6 z m 
m 

< z 


(B.6a) 

(B . 6b) 
(B.6c) 


The radial electrical field is given as (Frost et al . 1981) 
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E(r) = 


(B.7) 


-^r- r when r < R 
2 o 

R2p e - 

o ~- r when r > R 

2re o 

where R is the radius where the radial electrical field is a maximum, 
E (z). Equation B.7 can be rewritten as 


E(z,r) = 


^ E m (z ) when r < R 


Ejz) 


m 


when r > R 


(B.8) 


R can be approximately expressed as a function of z by the relationship 
R(z) = R. + 0.333 z (B .9) 

J 


where R. is the nozzle radius. Equation B.9 is estimated by interpo- 

U 

lating the edges of the voltage field shown in Figure 3.9. The value 
of E m is not constant in the jet field. The variation of E m with height 
is shown in Figure B.4. Again, these data are based on Ruhnke's analysis. 
A two region curve fit gives 


(10 ( z / z m ) ; z < 0.1 z 

1 rmax m - m 


E ( z ) = 
m 


rmax 


l 1 + f < z ' z m - °' 1 ) 


; o.i z m < z 


m 


(B.10) 


where E rmgx is the maximum radial electrical field strength. 

The expressions for electrical field strength (i.e., Equations 
B.10, B.8, B.6, and B.5) and for the wind field (i.e.. Equations B.4, 
B.3, and B.2) are substituted into Equation B.l. Equation B.l can be 
solved numerically, and the particle trajectories can be calculated 
based on given particle radius, jet exit velocity, U Q , and electrical 
field strength, E(0), and E rmax - 
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The input variables for solving Equation B.l are based on the 

reported measured data. The electrical field strength at the nozzle 

exit E(0) is equal to 10^ v/m. And the maximum radial electrical field 

strenqth, E , is estimated as 10 v/m. A uniform vertical electrical 
3 rmax 

field of E = 100 v/m, which represents a fair-weather background, is 
assumed to exist. From Appendix C, the particle diameter d is estimated 
to be in the range from 0.4 to 0.8 y. A value of d - 0.8 y is used. 

The number of charges on each particle is equal to 327 e" (this is the 
saturation charge). A nozzle exit velocity U Q = 300 m/s is also used. 

The computed particle trajectories based on the above data (which is 
representative of the current prototype charged particle generator 
system) are shown in Figure B.5. The origin of the particle trajectories 
are at the exit of the jet and at 0.001 m increments across the nozzle 
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Height (m) 




exit, respectively. The figure shows that the 0.8 y particle with 
saturation charge spreads outward to a distance on the order of 0.65 m 
and reaches heights of 80 m. Figure B.6 illustrates the computed 
particle trajectories in a free jet without an electrical field. Com- 
paring Figures B.5 and B.6, the charged particles are seen to spread due 
to their electrical charge but not to exit the jet and return to ground 
as often hypothesized in the literature. 

Realistically, the particle may not reach heights of 80 m because 
both the meteorological crosswind and the turbulent fluctuation in the 
jet are neglected. One can conclude qualitatively, however, that if the 
momentum of the jet flow can overcome the opposing electrical field 
( i . e . , below z m in Figure B.2) and thus push particles beyond z m , the 
charged particle will then rise due to their induced field until the 
gravitational force becomes greater than the electrical force and the 
jet drag force. 
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Figure B.6 Particle trajectories in a 
field. 
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APPENDIX C 


ESTIMATE OF WATER DROPLET SIZE 


The charge- to-mass ratio based on measured data is 

C 20 x TO -6 amp ^ 0 f n i\ 

m = ^ = 0.2 coul/kg (C.l) 

(6 cc/min)(1/60 min/sec)(10 P kg/cc) 

The Rayleigh limit representing the maximum number of charges a droplet 
can acquire is given by 

qp $ = 2.08 x 10 9 r 2 E (for water droplets at 20° C) (C.2) 

This represents an upper limit on the amount of charge a droplet can 

obtain. From Equations C.l and C.2 the radius can be computed for a 

given value of E . 

3 o 

The value of E Q , estimated from dividing the measured voltage drop 
across the corona by the separation distance between the needle and 

r 

attractor, is roughly 10 v/m. Hence, the number of droplets per 
kilogram of water injected is 


Number of Particle 
Kilogram of Mass 



= #/kg 


(C .3) 


The number of charges per particle becomes 

Number of Charge £ f4_ 3 

Particle M [3 71 p j 


(C .4) 


Equation C.4 must be less than or equal to the Rayleigh limit (i.e.. 
Equation C.2). Hence, 


(0.2 c/kg) 


6.25 x 10 18 

c 


irr 3 x 10 3 
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< 2.08 x TO 9 x r 2 x 10 6 



The particle radius is therefore predicted to be on the order of r = 

0.4 x 10" 6 m or the diameter is d = 0.8 y . 

Since the droplet is only in the corona region a very short time, 
it may not be realistic to estimate the size of water droplets based on 
the Rayleigh limit. According to Crawford (1976), the rate of droplet 
charging in the corona is approximately 


q P _ t 

V *0 + * 


(C .5) 


where q is charges on the droplet and q is saturation charge on the 

r r ^ 

particle is given by Equation C.2. t is the time constant of this 
process defined as the value of t when q p = 0.5 q . t Q is given by 


4e 

t - 2 

o Nq k 


(C.6) 


where N is ions density in the corona region, k is the mobility of the 
ions (k = 2.2 x 10“ 4 m 2 /v s). 

The charging time t is estimated by dividing the nozzle length 
with the airflow (essentially the droplet speed) velocity V, i . e . r 

t = ~ (C .7) 

Introducing measured data into Equation C.7, t is found to be approxi- 
mately 5.67 x 10" 4 sec (V = 300 m/s is assumed). From Equation C.5 

Vv ■ °- 53 

15 3 

The ion density used in Equation C.6 was taken as N = 2 x 10 ions/m 

(see Crawford (1976)). qp S is given by Equation C.2. Introducing q^ = 

0.53 q , into Equation C.4, the droplet radius is found to be on the 
ps 

order of r = 0.2 p or the diameter is d = 0.4 y. 
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